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Abstract 



One of the main challenges of the Tevatron at Fermilab and the Large 
Hadron Collider (LHC) at CERN is the determination of Standard 
Model (SM) parameters at the TeV scale. In this context various 
processes will be investigated which involve bottom-quark jets in the 
final state, for example decays of top-quarks or gauge bosons of the 
weak interaction. Hence the theoretical understanding of processes 
with bottom-quark jets is necessary. In this paper analytic results will 
be presented for the weak corrections to bottom-quark jet production 
— neglecting purely photonic corrections. The results will be used to 
study differential distributions, where sizeable effects are observed. 



I. Introduction 



With the start of the LHC a new energy regime is accessible, either to confirm the Stan- 
dard Model (SM) or to verify new physics at the TeV scale. Famous possible SM ex- 
tensions are heavy gauge bosons (e.g. Z'), Supersymmetry or Kaluza-Klein resonances 
in models with extra dimensions. Beside the outstanding discovery of the Higgs bo- 
son, processes involving top-quarks, gauge bosons of the weak interaction and jets are 
of particular interest. The experimental identification will rely on their characteristic 
decay products with leptons or bottom jets as characteristic examples. In particular 
in processes including top-quarks or a (light) Higgs boson the identification of the 
bottom-quark jets will play a crucial role. In addition to these SM processes, bottom- 
quark jets are also involved in many decays originating from signals for physics beyond 
the SM. In particular new resonances, e.g. heavy gauge bosons, decaying into ^-jets 
require a detailed SM-based prediction to prove a deviation from the SM. The discrim- 
ination of Z?-jets from light quark- and gluon-jets, ^-tagging, makes use of the large 
lifetime of B-mesons. In the last years there were large efforts by several experimental 
groups to develop new or to improve established ^-tagging algorithms. This requires 
a detailed theoretical understanding of the corresponding processes like bottom-quark 
pair and single bottom-quark production. Both processes were studied in the past. The 
differential cross section for bottom-quark pair production is known to next-to-leading 
order (NLO) accuracy in quantum chromodynamics (QCD) 03 EJ. For massless single 
bottom production the NLO QCD corrections can be extracted from di-jet production 

mm. 

It is well known that weak corrections can also be significant because of the presence of 
possible large Sudakov logarithms. These effects were studied intensively for several 
processes like weak boson, jet and top-quark pair production BH [6l Ul IS1 13 [10l [TQ El 
[T3l . Earlier work on Sudakov logarithms in four-fermion processes can be found in 
Refs. lUSIElIiainillSllElESlEIllia.A first study for bb production can be found 
in Ref. |l24l . where the partonic sub-processes gg — > bb and qq — > bb were considered. 

In experimental studies, e.g. of top-quark pair production at the Tevatron, often only 
one Z?-jet is required. Considering the importance of these background studies also 
single bottom-quark production will be studied in this paper. 

The outline of the paper is as follows: In Section |II]the leading-order cross sections for 
the production of one and two bottom jets are presented. The partonic contributions 
are split into quark-, gluon- and pure bottom-induced processes and relations based 
on crossing symmetries are introduced. After a discussion at the hadronic level we 
find, that leading order processes involving electroweak gauge boson exchange can 
be neglected. Concerning the differential pj distribution at high energies the effects 
from pure bottom-induced contributions are also negligible. In Section [Hi] we present 
the virtual and real corrections of order a^a for the remaining processes. The virtual 
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contributions to quark-induced processes contain infrared and ultraviolet singularities, 
while the gluon-induced ones are infrared finite. Compact analytic expressions for the 
various channels are presented. In Section [IV] we present various consistency tests of 
our calculation and discuss the numerical results at the hadronic level. Our conclusions 
are given in Section IVl 



II. Bottom-jet production at leading-order 



At parton level three types of processes will be distinguished: Quark-induced pro- 
cesses (with two quarks in the initial state, one of these being u, d, c, s), gluon-induced 
contributions (with one or two gluons in the initial state) and finally pure bottom-quark 
induced processes: 

quark — induced : qq — > bb, bq — > bq, bq — > bq, bq — > bq, bq — > bq, 

gluon — induced : gg — > bb, bg — > bg, bg — > bg, 

bottom — induced : bb — > bb, bb — > bb, bb — > bb. (II. 1) 

Processes with initial state photons are neglected. Here the light quarks (denoted 
generically by q and q) and the bottom-quark (b and b) are taken as massless. Sample 
diagrams are shown in Fig. III. II and Fig. 111.21 The colour and spin averaged squared 
matrix element Y] M l 'j^ k ' l \ 2 can be written either in terms of the Mandelstam variables 
s, t and u or as a function of s and the cosine of the scattering angle cos f> = z = (w — t)/s 



i,j^k,l\2 



(s,t,u) =Y}M 



ij->k,l\2 



(11.2) 



Various explicitly checked crossing relations will be useful below. For the quark- 
induced processes in Eq. ( III. II) we find 
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and similarly for gluon-induced and pure bottom-induced processes 
T\<M bb ^ bb \ 2 



Y\M bb ^ bb \ 2 = Y\M bb ^ bb \ 2 

The spin and colour averaged differential partonic cross section reads 

do i,j->k,l 2 2 



dz 



1 32ns 



(H.3) 

(H.4) 
(H.5) 

(H.6) 
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where 1 = 2 for identical particles in the final state and 1=1 otherwise. 
Before presenting explicit results we complete the notation and definitions used in this 
and the following Sections. As usual a s and a stand for the strong and electromag- 
netic coupling constant. The Z boson mass is denoted by mz and the (axial-)vector 
coupling is expressed in terms of the sine (cosine) of the weak mixing-angle s\v (cw), 
the weak isospin Tj and the electric charge Qf for a fermion of flavour / in units of 
the elementary charge e: 

<T 3 f -2s w 2 Q f ), g f a = ^—T{. (II.7) 



2s\yCw 2sy/Cw 

The Cabibbo-Kobayashi-Maskawa mixing matrix is set to 1 and the W coupling to the 
quarks is 

Sw = (H.8) 
2\/2sw 

Let us start with the processes listed in Eq. (III. II) , which involve only quarks in the 
initial state. These can be mediated by gluon, Z boson and photon exchange. Conse- 
quently we seperate pure QCD 0(a 2 ), electroweak 0(a 2 ), and mixed QCD-electro- 
weak 0(a 5 a) contributions. In the case of one or two gluons in the initial-state an 
interaction via electroweak gauge bosons is not possible at leading order and only 
QCD processes contribute. We start with the QCD-contribution for quark-antiquark 
annihilation into a bb-pak 

£|^r^| 2 (^( M , M )) =4ay^l(l+z 2 ), (n.9) 
with N being the number of colors. The purely electroweak contributions are given by 

El<r1 2 (^(^,«)) = i6aV[eie 2 (i+z 2 ) 

„2 



+ {s _ mz 2 ) 2 {(Sr + gf) (g° V + fa ) ( 1 + ~ ZgWagaZ 
S — M7 V 



(11.10) 



Because of colour conservation the interference between the QCD and electroweak 
matrix elements vanish. The remaining quark-induced contributions can be obtained 
using the crossing relations in Eq. (111.31) . 

The gluon-induced processes in Eq. (III. II) proceed through QCD amplitudes only and 
the corresponding crossing relations are defined in Eq. (111.41) . For the gluon fusion 
channel the squared matrix element is 
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Finally we consider the pure bottom-quark induced contributions arising from pure 
QCD, mixed QCD-electroweak and electroweak matrix elements. Below we list the 
results for bb — > bb scattering, those for bb — ► bb and bb — > bb scattering can be de- 
duced from Eq. (111.51) . 

,N 2 -1 



bb^bb\2 

■2 



s,z(s,t,u) 



4a 2 7t 2 
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(11.12) 



_5 (l-z) 2 2s(l+z)+m z 2 (l-z) 

-m z 2 l+z s(l +z) +2m z 2 



(11.13) 



' 4 iV( 1 1 - 2z + 4z 2 + 2z 3 + z 4 ) - 2 ( 1 - z) ( 1 - z 2 ) 
Ub N{l+z) 2 



+ 



N(s-m z 2 ) 2 (s(l+z)+2m z 2 ) 2 
x \gf + g b a 4 ) (2m z 4 (2(l -z) 2 + ^V(7 - 2z + 3z 2 )) 

- 2sm z 2 (l-z)((l-z) 2 + 2JV(4 + z + z 2 )) 

- 5 2 (2(1 +z)(l -z) 2 -W(ll -2z + 4z 2 + 2z 3 +z 4 )) 

+ gf gf (4m/(6(l -z) 2 + /V(l - 14z + 5z 2 )) 
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+ 25(l+z)((l-z) 2 +M3 + z 2 )) 



(11.14) 
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The term proportional to aa s originates from the interference between s and t exchange 
amplitudes and is present for the process bb — > bb and its crossed versions only. The 
purely electroweak contributions proportional a 2 are found to be negligible. 
With these ingredients we calculate the hadronic transverse momentum (pj) distribu- 




Figure II. 1: Sample Born diagrams for quark- induced bottom-quark production. 




Figure II.2: Sample Born diagrams for gluon-induced bottom-quark production. 

tions for single and double b-tag events and investigate the relative importance of the 
different processes described above. As a consequence of the massless approximation, 
most partonic contributions are ill-defined in the limit z — > ±1, where initial and final 
state parton become collinear. Since Z?-jets close to the beam pipe escape detection, 
we require a minimal transverse momentum (/?j ut ) of 50 GeV. The leading order cross 
section is obtained from 
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do Hu H 2 ^x = Y,f 6x1 J dx 2/f/H 1 (xi)/ 7 y^ 2 (x 2 ) 

x d6i^ x (xiP h x 2 P2) ®(pr>pf t ), (H.15) 

where the factorization scale dependence is suppressed and xi and X2 are the partonic 
momentum fractions. The parton distribution functions (PDF's) for parton i in hadron 
H are denoted by fu H . The sum runs over all possible parton configurations in 
the initial state. 

For the study of Z?-jets we will distinguish between single £>-tag and double b-tag 
events, and evaluate the pj distributions with the following input parameters 

a, = 0.1, a = -!-, m z = 91.1876 GeV, m w = 80.425 GeV. (11.16) 
128 

For sw we use the on-shell relation. Moreover we used the PDF's from CTEQ6L with 
the factorisation scale set to /uf = 2p™ 1 . 




Figure ET.3: Differential cross section as a function of pj for single &-tag events (upper 
left figure) and double b-tag events (upper right figure) at the Tevatron (y/s =1.96 
TeV) and the relative composition normalized to the Born cross section (lower figures) 



Let us first discuss the results for the Tevatron (y/s = 1.96 TeV). The composition 
of the corresponding leading-order differential cross section for single b-tag events is 
shown in Fig. III.3f a). In the upper plot the cross section is depicted as a function of 
the transverse momentum pj, its relative composition in percent is shown below. For 
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Figure II.4: Differential cross section as a function of pj for single b-tag events (upper 
left figure), double b-tag events (upper right figure) at the LHC (y/s = 14 TeV) and the 
relative composition normalized to the Born cross section (lower figures) 



px-values up to 200 GeV the gluon- and quark-induced QCD processes are of the same 
order, while for higher energies the quark-induced b-)et production is dominating the 
px-distribution. Contributions from gluon-induced processes drop fast with increasing 
pj and contribute only a few percent to the differential cross section for pi > 300 
GeV. The relative contributions from bottom-induced processes and processes with 
electroweak boson exchange are only of the order of few permille and therefore negli- 
gible for the px-values considered. 

For pr < 100 GeV and in the double b-tag case the px-distribution (Fig. HI .3 l b)) is 
dominated by the gluon fusion channel. For higher pj-va\ues QCD quark-antiquark 
annihilation takes over and for transverse momenta larger than 200 GeV it dominates 
completely. The relative contributions from processes involving electroweak bosons 
are of the order of 1% for pj < 80 GeV and, for pj > 260 GeV are of the same order 
(few permille) as the gluon-fusion process. Bottom-induced contributions yield a few 
percent to the differential cross section for pj < 100 GeV and become insignificant for 
higher transverse momenta. 

For the LHC the results are significantly different. The LO /?x-distribution for sin- 
gle b-tag events is presented in Fig. HI4f a). The gluon-induced processes dominate 
in the low energy regime (px < 500 GeV). For pi larger than 500 GeV the quark- 
induced processes take over and finally dominate the distribution. This "cross over" of 
gluon- and quark-induced contributions is a consequence of the different behaviour of 



7 



LHC quark and gluon luminosities. The relative contributions from the exchange of 
electroweak bosons and bottom-induced processes are always below 1% and therefore 
negligible. A similiar picture is observed for the differential cross section for double b- 
tag events at the LHC (Fig. III.4r b)). Here the "cross over" of quark- and gluon-induced 
contributions is around pj = 1 TeV. For low pj the pure bottom-induced processes are 
as important as quark-induced contributions. This seems surprising, because the par- 
ton luminosities of bottom-quarks in a proton should be highly suppressed compared 
to the light flavours. There are two reasons for the relatively large cross section of the 
purely bottom-induced processes. First, the partonic cross sections of bb and bb scat- 
tering are strongly enhanced for large z values. In this region the parton processes with 
bottom-quarks in the initial state can be several orders of magnitude larger than the 
quark-antiquark-induced process. Second, the bottom-quark PDF is essentially ob- 
tained by multiplying the gluon distribution in the proton with the splitting function of 
a gluon into a bottom-quark pair. Because of the high gluon luminosities at low ener- 
gies, the bottom-quark PDF becomes of the order of a few percent relative to the PDF's 
of the light flavours. This, together with the large partonic contributions is responsible 
for the relatively large bottom-induced differential cross section. The argumentation 
implies that the main contribution from bb, bb and bb scattering comes from the low 
pi region, while for high pj values these effects are small. It might be interesting to 
study whether the b-PDF could be further constrained using bb production at low pj. 

As shown above, the leading order contributions from electroweak gauge boson ex- 
change are always negligible for the study of Z?-jet production. This is also true in 
the context of NLO corrections with an expected size of serveral percent relative to 
the leading order distributions. Moreover, we have shown that the QCD contributions 
from processes with two bottom-quarks in the initial state are unimportant for the study 
of px-distributions at large transverse momentum. In particular with regard to the Su- 
dakov logarithms becoming important at high energies this approximation is justified. 
Consequently the weak 0(a) corrections to bb — > bb, bb — > bb and bb — ► bb will not 
be discussed in this article. 
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III. Weak corrections to bottom-quark production 



In this Section we calculate the weak corrections of order oc^oc to the following partonic 
processes 

qq — > bb, bq — > bq, bq — > bq, bq — > bq, bq — > bq, 
gg -> bb, bg -> bg, bg -> bg, 

neglecting photonic corrections. We subdivide the 0(a) corrections in contributions 
from quark-induced processes (Section IIII.ll) and gluon-induced processes (Section 
IIII. 21) . Before presenting analytic results for the weak corrections, let us add some 
technical remarks. For the calculation of the next-to-leading order weak corrections 
the 't Hooft-Feynman gauge with gauge parameters set to 1 is used. The longitudi- 
nal degrees of freedom of the massive gauge bosons Z and W are thus represented by 
the goldstone fields % and (]). As mentioned in Section HH all incoming and outgoing 
partons are massless and consequently there are no contributions from the Goldstone 
boson % and the Higgs boson. Ghost fields do not contribute at the order under consid- 
eration. For the analytic reduction of the tensor integrals to scalar integrals we used the 
Passarino-Veltman reduction scheme 11251 . The scalar integrals were calculated either 
analytically with standard techniques or evaluated numerically using the FF-library 
[26]. The convention for the scalar integrals in this article is: 

in 2 J (£ 2 -m 2 + ie)--- 



The bare Lagrangian L is rewritten in terms of renormalised fields and couplings as 
follows 

1/2 1/2 

L (^0,^0,^0,^0) = ^{Z4 ^R> Z A A R ,Z m m R ,Z g g R ) 

= L( x ¥ Rl A Rl m Rl g R )+L ct ( x i' Rl A R ,m R ,g R ) (111.2) 

(see e.g. Ref. Il27ll ). For the present calculation only wave function renormalization 
is needed and no mass or coupling renormalization has to be performed. The partonic 
NLO corrections are, furthermore, independent of the factorization scale /Uf. The wave 
function renormalization is performed in the on- shell scheme 

VjAL = ^R ) L)1/2^R J L = ^ 1+ 1 5Z R,L^R ) L (ffl 3) 

where the renomalization constants are given in terms of self-energy corrections E and 
their derivatives: 

8Z V = ^(5Z L + SZ R ), 

5Z A = ^(5Z L -5Z R ). (III.4) 
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Here only SZy is required. For bottom-quarks it is given by ETll : 
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(III.5) 



corresponding to the massless limit of the formulae given in Ref. ffT2l after the ex- 
change of bottom and top mass. For light quark flavours the renormalisation constant 
reads: 
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(III.6) 



In the following, most results will be given using the explicit formulae of the scalar 
integrals. To evaluate the NLO corrections to the differential cross section only the real 
part of the virtual corrections contribute. Therefore we will skip the parameter for an- 
alytic continuation £ wherever possible. Hence the relations from crossing symmetries 
lead for example to the following replacements (s > 0, t < 0, M 2 > 0) 




III.l. Weak corrections to quark-induced processes 

We start with the 0(a) corrections to quark-induced processes, consisting of vertex-, 
box- and real-contributions. In contrast to the gluon-induced reaction where elec- 
troweak corrections to the QCD Born amplitude can be clearly identified, the cor- 
responding classification for qq — > bb more precisely of all contributions of 0(a 2 a) 
is more involved. In this case QCD box amplitudes 0(a,, 2 ) may interfere with the 
weak amplitude of 0(a) and similarly, mixed box amplitudes of 0(a s a) may interfere 
with the QCD Born amplitude of 0(a s ). Infrared singularities are cancelled by con- 
tributions from interference between initial state radiation and final state radiation. A 
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closely related discussion can be found in Ref. ifTTI and for QED corrections to neutral 
current corrections in Refs. Il28ll29ll . Two types of box-diagrams can be distinguished: 



1. The (box-type) weak correction to the QCD Born amplitude, interfering with the 
QCD Born amplitude. 

2. The QCD box diagram interfering with the weak Born amplitude. 



Sample diagrams are shown in Fig. IIII.ll b). c). The box-diagrams are UV finite. Their 
IR singularities cancel against those from real emission, specifically, from interference 
terms between initial and final state radiation. The vertex-corrections are infrared (IR) 
finite and their UV divergencies cancel against the wave function renormalization de- 
scribed above. All contributions are free from initial state mass singularities. To handle 
the infrared singularities we use the dipole subtraction method 113011 . In the notation of 
Ref. 1301 the NLO box and real corrections can be written as 



g a NLO O(a^a) 



□+Real 



dc D + (l ® do B ) 
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+ 



da 
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52 dVDipole 
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da* 



J <1x(k(jc)+P(jc)) ® do B (x) 



(HI.7) 



Here da D consists of the virtual box corrections. Their IR- singularities will be can- 
celled by the contribution involving the I-operator. Note that the differential cross 
section da B is not the QCD leading order contribution. In order to respect the mixed 
QCD-weak structure of the box-diagrams and of the real corrections, da B is given 
by the interference between the QCD amplitude times the weak amplitude, where the 
colour correlations of the dipole formalism yield a non-vanishing colour factor. We 
use the symbol '£§>' in Eq. (IIII.7D to denote color as well as spin correlations. The K 
and P operators are finite remainders originating from dipoles, for details we refer to 
Ref. AH. 

For the virtual and real corrections we present analytical results for the quark-antiquark 
annihilation process, while for the remaining quark-induced processes we refer to the 
crossing relations in Eq. (111.31) . The results for 

(I + K + P)®da B (III. 8) 

will be given explicitly for each process. The virtual contributions da Virtual to the 
partonic differential cross section is given with Eq. (|II6I) 
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The initial- state vertex corrections for qq — > bb read 
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(III.ll) 



and the squared leading order matrix element is given in Eq. (111.91) . The final state 
vertex corrections read 
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where the scalar integrals are defined in the Appendix [A] For the box-diagrams we 
find 
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The factor B q v~* bb results from the interference between the leading order QCD am- 
plitude and the weak amplitude. Terms proportional to the Z vector couplings are odd 
functions in the scattering angle z as a consequence of Furry's theorem. The virtual 
corrections for the remaining quark-induced processes can be deduced by the cross- 
ing relations Eq. (III. 3 1) . For latter use we extract the IR divergent contribution to the 
differential cross section 



da 



□ IR 



qq—hh _ OCa 2 N 2 — 1 ^qq^bb 



32s N 2 



(4tc)T(1+£) 



-ln(-)+ln z ( — 



hT[ ~t 



(III. 18) 



The contribution from real emission is given by 



qq->bbg 



a 2 a(4rc) 



N 2 -l 
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(g q v g b v {ti + 1\ + u\ + u\) - g q a g b a {t 2 + 1\ -uj- u\) 



1 1 



1 



1 



S S-m Z 2 S + t\ +?2 + Wl +U2 S + ti+t2 + Ui+U2+mz 2 

1111 



S + tl+Ul S + t2 + U{ S + tl+U2 S + t2 + U2 

x (2s 2 + (t l +t 2 + u 1 +u 2 )(2s-mz 2 )^ 
x +t 2 - u\ - u 2 )s 2 + {{t\ + 1 2 ) 2 - {u\ + u 2 ) 2 )s 
+ (H+t2 + U\+U2)(ht2-UiU2) 



(III. 19) 



with 



ill 



Ikqkjj, ?2 — ^.kqkfo, 
2-kqkjj, U2 = Ikqkjj. 



(111.20) 



The combination s + 1\ + 12 + u \ + 112 + niz 2 appearing in Eq. (111119b can vanish, cor- 
responding to on-shell production of the Z boson. Therefore we introduced according 
to the principal value prescription a cut around the singular point 



8> \ \s + t\ +t 2 + ui +U2 + m z 2 \. 



(111.21) 



and demonstrated numerically that the result remains unchanged for 8 between 10 -1 
and 10~ 5 GeV. Moreover we checked that also a fictitious Z boson width 0.01 GeV < 
Yz < 1 GeV produces the same numerical result. In analogy with the leading-order 
crossing relations we find for the remaining real corrections: 



qb^qbg 



qb->qbg 



qq-^bbg 



(111.22) 



S+-*tl, t 2 ^-S-t 1 -t 2 ~U l -U2 



qb^qbg 



E 




^qbg 


E 




^qbg 


E 




->bbg 



S— H\, t\^>U\, t 2 -^>U2, U\-^rS, U2-^—s—t\—t2—ui—U2 

(m.23) 



For these contributions the Z boson remains off-shell. To get an infrared finite result 
the corresponding subtraction terms from the dipole formalism were implemented. 
We checked explicitly the numerical stability of real corrections and dipoles and found 
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the pointwise cancellation between the two contributions in the singular phase space 
regions. However the dipoles approximate the real corrections very well also in non- 
singular regions. This results in large cancellations in the integrated result. In contrast 
to the singular configurations this cancellation is not pointwise but takes place between 
different phase space regions. As a consequence of this behaviour one needs high 
statistics for the numerical integration. To check the numerical results we therefore 
implemented also a variant of the phase-space-slicing method [2] and find agreement 
comparing the results of both methods. The relevant formulae for the slicing method 
and the plot for the comparison are shown in the Appendix El 



Let us now come to the contributions from the dipole formalism F ® da B with F = 
I, K,P. The symbol '®' denotes spin and colour correlation between the operator F 
and the leading order amplitudes, which are treated in general as vectors in colour and 
spin space. For the processes under consideration, the gluon is always emitted from a 
fermion line and only trivial spin correlation appears. We start with the contribution 
for the quark-antiquark annihilation channel 



(l®do B ) 



oca 2 N 2 - 1 
' 32* N 2 

,2 



+ InM^ 

-u 



*bb 



(47!) 



r(i 



e \u 



In — +31n - 



(111.24) 



where B qq ^ hb is defined in Eq. (IIII.17I) . This contribution to the differential cross 
section cancells the IR-poles and the /j-dependence from the virtual corrections in 
Eq. din. 181) . For the finite K and P parts we find 



(P®do B )^(*) 



0, 

oca, 2 N 2 - 1 
32sx N 2 



B qq ~ 



mil 

u 



l+x z 
\-x 



(111.25) 
(111.26) 



Of course in Eq. dIII.261) B qq ~^ bb is also a function of x and for the evaluation of the 
plus distribution the well known relation 



f Q dx(f(x) 



+ 



dxf(x)[g(x)-g(l)\ 



(111.27) 



is used. For bq — > bq and bq — ► bq we find 
(l®do B ) hq ^ bq - 



(l®da B ) 
aa? A^ 2 — 1 



32^ A^ 2 



ft^bq->bq 



(47C) 



r(l-e) 



■In 



s 

—u 



16 



+ { +31*' s 



— u 



— u 



(111.28) 



(K®do B ) bq ^ hq (x) = (K®da B )"(x 

aa? Af 2 - 1 



+ 



32sx N 2 
3 / 1 



B ^bq^bq 



8(1-,) [ \-± 



2 V 1 -x 



, pln(l-x) 
\ 1 —x 



(1+x) ln(l-x) 
(111.29) 



(P®da B )^(x) = (P®do B )"(x 



B *.bq^bq 



aa 2 N 2 



32sx N 2 



-B^^ln (-^) (l^r) ' (IIL30) 



where B bq ^ hq can be deduced from B qq ^ hb and Eq. (111.31) . For &g — > bq (and &g — > bq) 
we find 



(l®da B ) 



■Q^bq^bq 



(l®do B ) 

aa 2 N 2 - 1 
325 N 2 

,2 



B^bq^bq 
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(111.31) 



(K®da B )" rt (x) = (K®da B )"W 

aa? N 2 - 1 



32s* N 2 



R\bq->bq 

-gbq^bq 



«(l-x)(|4 



+ 



3 / 1 \ /21n(l -x 



+ 



2 \\ —x J , \ \ —x 



+ J 



(1+x) ln(l-jc) 

(in.32) 



(P Ode 8 ) = (P®da B )"W 



K^bq^bq 



aa 2 N 2 - { B bq^bq ln (z^ ( 1 +* 2 



32sx N 2 



l-x 



(111.33) 



+ 



and B bq ~^ bq is obtained by the crossing relations Eq. (111.31) as above. 
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(a) 



(b) 
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Z,W,(|) 
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Figure EH.2: Sample diagrams for the gluon fusion process at 0(a 2 a). 



III.2. Weak corrections to gluon-induced processes 



Sample diagrams for the virtual corrections to the gluon-induced processes are shown 
in Fig lffl.21 No infrared divergencies are present and no corrections from real radia- 
tion contribute. The weak corrections are again classified into those from self-energy-, 
vertex- and box-diagrams. The latter are UV finite, vertex and self-energy contribu- 
tions must be renormalized. The differential cross section for the gluon fusion channel 
at next-to-leading order is decomposed as follows: 



dz 



1 



32715 



7=Z,W,( 



Vi{gg->bb)\2 



(111.34) 



We start with the analytic results for the self-energy diagrams (Fig IIII.2r c)) in the 
gg — > bb channel, which can be written in a compact form 



£|^£z(gg-/^)|2 = %a 



2 , 2-N 2 +N l Z, b 2 , b 2, l+z 



N(N 2 -\)s 



s (l+z) 
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(111.35) 
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(111.36) 
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where 



fiyuyi) = -yiin ( j^jj -y 2 in ) (m.37) 



t + m w 2 -m t 2 ( I 4miv 2 ? 
yi/2 = x|l±Wl- 



2t \ V {t + m w 2 - m t 2 ) 2 



(111.38) 



The contributions involving (|) and are closely related 

£|^(gg-^)|2 = _^!_£|^^fe^)|2_ (m39) 

The vertex corrections consist of 5, ? and u channel contributions. Since the corrections 
to the s channel are proportional to the quark mass (see Ref. [12J Eqs.(II. 18-22)) they 
vanish in the massless limit. The remaining vertex corrections (e.g. Fig. |III.2r b)) can 
be expressed in terms of the self energies 

£|^,.(gg-M>)|2 = _ 2 £|5tf E <'(s£^)|2. (IH.40) 

The box contributions (Fig IIII.2r a)) are more involved and are listed in Appendix O 
Corrections to gb — ► gb and gb — > gb are obtained via the relations Eq. (|II.4I) . 



IV. Results 

Apart of the obvious checks, like cancelation of UV and IR singularities, a numerical 
test was performed for the IR-divergent contributions, implementing the phase- space- 
slicing method. As illustrated in Appendix B we found complete agreement of the 
two methods. For the quark-antiquark annihilation and the gluon fusion channel we 
found complete analytical agreement with previous results for top-quark production 
lfT2l . with the following replacements 

m t ^m b , g' v ^g b v , g'a^Sa and m t ^0. (IV. 1) 

These results also allow for a numerical comparison between the massive and the mass- 
less approach. In Fig. IIV. li the results for a massive bottom-quark (rn\, = 4.82 GeV) are 
compared with those from the massless approximation. Agreement between the mas- 
sive and the massless results is observed in the full kinematic range, which justifies 
the massless approximation for the bottom-quarks used in the paper. The difference of 
the massless and the massive result concerning the relative corrections is always below 
0.5%. Furthermore all LO and NLO partonic corrections were calculated analytically 
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and thus all crossing relations listed in the previous Sections are used as additional 
cross check. For the numerical evaluation we use the same input parameters as in Sec- 
tion HIl As mentioned in Section [III] the results are leading order in QCD and therefore 
evaluated with the leading order PDF's CTEQ6L. 

We start with the results for the differential px-distribution at the Tevatron. In Fig. HV.2l 
the relative corrections normalised to the full leading order cross section are shown. 

Double b-tag 





-5 



-10 



-15 

500 1000 1500 2000 
P T [GeV] 

Figure IV. 1: Relative weak corrections for double b-tag events at the LHC, comparing 
the results for a massless bottom-quark with the massive case. 

For single b-tag events (upper figure) the relative corrections are negative and of the 
order of a few permille up to pj = 120 GeV. The following increase of the relative 
corrections with the resonant behaviour around pi = 180 GeV is caused by the vertex 
corrections to qq — > bb involving virtual top-quarks. In this range of pj the weak cor- 
rections become positive and increase the differential cross section by a few permille. 
Above the resonance peak the weak corrections become negative again and the typical 
Sudakov behaviour is observed. The magnitude of the relative corrections increases 
with energy for pj > 200 GeV and reaches —3% for pj = 400 GeV. 
Due to the aforementioned virtual top-quarks, the relative corrections for double b-tag 
events are positive for pj values below 220 GeV (lower figure). In this case qq — > bb 
and gg — ► bb are the only contributing partonic channels, hence the resonance peak is 
more distinct than in the single b-tag case with the maximum (+1.5%) of the relative 
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Single b-tag 
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Double b-tag 
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Figure IV. 2: Relative weak corrections for single b-tag (upper figure) and double b-tag 
events (lower figure) at the Tevatron. 
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Single b-tag 




100 200 300 400 

p T cut [GeV] 




Figure IV.3: Relative corrections to the cross section for pj > p^ ut at the Tevatron for 
single £>-tag (upper figure) and double £>-tag events (lower figure). The red lines give 
an estimate on the statistical uncertainty as described in the text. 
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corrections around pj = 1 80 GeV. The influence of the Sudakov logarithms starts for 
pj > 220 GeV, where the weak corrections change sign and become negative. For 
pj = 400 GeV the relative corrections amount up to —3%. 

The relative corrections to the integrated cross section (px > Pj 11 ) at the Tevatron are 
shown in Fig. IIV3[ For the single b-tag (upper figure) and the double b-tag (lower 
figure) the peak around the Jf-threshold is smoothed by contributions from higher pa- 
yables, leaving the integrated correction negative in the single b-tag case. Beyond 
pj = 200 GeV the relative corrections to the integrated distributions are between 
and —3%. In addition we give in Fig. IIV.3l a rough estimate of the expected statistical 
error at the Tevatron. The estimated number of events for pj > /? T ut is based on an 
integrated luminosity of 8 fb _1 . A comparison of the statistical uncertainty with the 
relative weak corrections shows that the weak effects will not be visible at the Tevatron. 

At the LHC however the corrections are significantly larger. For the subsequent disc- 
susion y/s =14 TeV will be adopted. The relative corrections to the leading order 
px-distributions are shown in Fig. IIV.41 For single b-tag events (upper figure) the rel- 
ative corrections are always negative and of the order of a few permille up to — 1% 
for 50 GeV < pi < 250 GeV. In contrast to the Tevatron results, the resonant be- 
haviour arising from virtual top-quarks in quark-anti-quark annihilation is not visible, 
a consequence of the dominance of the gluon-induced processes in this px-regime 
(Fig. |IL4r a)). For 250 GeV < pj < 1 TeV the weak NLO contributions vary between 
— 1% and —8%, compared to the leading order distribution, a consequence of the Su- 
dakov logarithms. In the high energy regime (px > 1 TeV) the relative corrections 
amount to —10% and will even reach —14% for pi = 2 TeV. 

The lower figure in Fig. IIV.41 shows the relative corrections for double b-tag events at 
the LHC. Despite the strong suppression of qq — > bb at leading order (see Fig. |II.4r b)) 
a small remnant of the enhancement from virtual top-quarks is visible in the double 
Z?-tag case. With increasing pj the Sudakov logarithms dominate the shape of the weak 
NLO contributions and yield relative corrections between —1% and —7% (250 GeV 
< pj < 1 TeV). At the highest px-values considered relative corrections up to —14% 
are observed. 

Fig. IIV.5I (upper figure) shows the integrated cross section for single b-tag events at 
the LHC together with an estimate of the statistical error based on an integrated lu- 
minosity of 200 fb _1 . The same composition in shape and magnitude is observed as 
for the differential distribution. The statistical error estimate matches the size of the 
weak corrections up to pj = 1.5 TeV. For higher pj-ya\ues the rate drops quickly and 
it will be difficult to observe the effect of the weak corrections. For double Z?-tag events 
(lower figure) we find again the smoothing of the ff-threshold for the relative weak cor- 
rections, while the composition of the curve for pj > 250 GeV is very similiar to the 
already discussed differential distribution. Considering the statistical error, the slight 
increase from virtual top-quarks will not be observable at the LHC. For px-values be- 
tween 250-1000 GeV the weak corrections are larger than the statistical error, above 
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Pt = 1 TeV they are comparable or smaller. 

In Fig. IIV.6l we show results for the LHC operating at y/s =10 TeV. The absolute cross 
sections are by more than a factor two lower, due to the lower parton luminosities. 
However the impact of the electroweak corrections is nearly the same and qualitatively 
similiar results are obtained for the relative NLO corrections. 



V. Conclusion 

In this article we present the weak corrections for &-jet production of order a^a ne- 
glecting purely photonic corrections. We derive compact analytic results for the NLO 
contributions in quark-anti-quark annihilation and the gluon-fusion. For the remain- 
ing partonic processes we list the crossing relations needed to derive the corresponding 
analytic results. We present the corrections to the px-distribution and compare the re- 
sults with the expected statistical uncertainty. In particular we find corrections up to 
ten percent for bottom-jets at high transverse momenta accessible at the LHC. These 
effects are larger than the anticipated statistical uncertainty. 
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Figure IV.4: Relative weak corrections for single b-tag (upper figure) and double b-tag 
(lower figure) events at the LHC. 
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Figure IV.5: Relative corrections to the cross section for pj > p™ 1 at the LHC for 
single £>-tag (upper figure) and double £>-tag events (lower figure). The red lines give 
an estimate on the statistical uncertainty as described in the text. 
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Figure IV.6: Relative corrections for single b-tag (a) and double b-tag events (b) at 
the LHC operating with a center-of-mass energy y/s = 10 TeV (upper figures) and the 
ratio of the differential leading order cross sections at the LHC operating at 10 TeV 
and 14 TeV (do LO (10TeV)/dp T ) / (da L0 (14TeV)/dp T ) (lower figures). 
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A. Notation 



We define as usual 



B (p 2 ,mi 2 ,m 2 2 ) = Jd d £- { 



(271//) 
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(£ 2 ~m 1 2 + ie)((£ + p i ) 2 -m 2 2 + iey 
Co(Pi,P 2 ,(pi+Pi) ,™i ,™2 ,m 3 ) = 



(271//) 



2e 
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For the UV-divergent one-point integrals Ao and two-point integrals Bo we define the 
finite parts through 

A (m 2 ) =ra 2 A + A ("- 2 



m 

B (p 2 ,m 2 ,m^) = A + B (p 2 ,m 2 ,m^), (A.l) 

with A = 1/e — Y+ In (An) . The 6-dimensional four point functions used in this article 
are 



D§(0,0,0,W,0,0,0,0) = ^ 2 (-}\ 

D^(0,0,0,5, M ,?,0,0,0,m z 2 ) = ~ — \(m z 2 -s) 
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Li 2 1 



+ In 



1 



In - 



+ 



s(t + m z 2 ) 



(A3) 
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The following abbreviations are used for the 4-dimensional four and three point func- 
tion 



Djf(z) = D (0,0,0, t,u,s,m z ,m t z ,m w z ,m t z ), 



(A.4) 



Co(s,m t z ,m w z ,m t 2 ) 
Co(s,m 2 ,m 2 ,m 2 ) 
Co(t,m t ,m t ,mw ) 



C (0,0,5,m f 2 ,miy 2 ,m ? 2 ), 
C (0,0,5,m f 2 ,m ? 2 ,m ? 2 ), 
Co (0, 0, t , m 2 , m 2 , m w 2 ) . 



(A.5) 



B. Phase space slicing 



In this appendix we list the relevant formulae for the phase-space slicing method ap- 
plied to the quark-induced processes listed in Section lTiLll The integral over the three- 
particle phase space can be written as 
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where we introduced a cut E c on the gluon energy E g . The second term in this equation 
is integrated in d = 4 dimensions, yielding a logarithmic dependence of the result on 
the cut E c . For the first term the real matrix element squared is calculated in the 
eikonal approximation and evaluated in d dimensions, using the parametrisation given 
in Ref. 0. For the eikonal factor for the qq — > bb process we find 
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(B.2) 



Adding the corresponding contribution from the virtual corrections Eq. (|M.18I) we find 



da s ? ft u - u 

qq-^>bb 
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(B.3) 
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2E 

with the dimensionless variable x m i n < For the remaining partonic processes the 
"soft" contributions to the differential cross section are given through 

a?a N 2 - 1 



do s ? ft , 
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do Soft 
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Finally we present the comparison between the phase space slicing and the dipol sub- 
traction method. In the figure below the difference for the differential cross sections 
at the LHC obtained with the two different methods in terms of standard deviations is 
shown. The agreement is always better than three sigma. 
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Figure B.l: Difference between results based on dipole formalism and phase space 
slicing for the differential cross sections for the different quark-induced channels in 
terms of standard deviations. For the comparison only the relevant IR contributions 
(boxes and real corrections) were taken into account. 
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Box diagrams for gluon fusion 
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2\2 



- —^—(5s 2 + As(3m w 2 -2m 2 )+S(m 2 -m w : 
1 — z V 

+ 25(5 — 2m, 2 + 2m l y 2 )z + 5 2 z 2 j 

x (Co(s,m 2 ,m 2 ,m 2 ) - (I + z)C (t , m 2 , m 2 , m w 2 )^ 

- 5(35 — 4m 2 + 4mw 2 ) 

1 — Z L 

+ (^25(25 - 3m f 2 + 4mvi/ 2 ) +6(m f 2 -mvi/ 2 ) 2 )z 

- 5 2 z 2 + 2^5(5-m f 2 + 2miy 2 ) + (m f 2 -m iy 2 ) 2 jz 3 
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x Co(s,m t ,m w ,m t ) 
1 



5s 3 + 22s 2 m w 2 - \Qs 2 m t 2 



2(1 -z) 

A A O O O O 1 

+ \6s(m t +2mw — 3m t mw ) — 16(m f — miy ) 
+ 5 ^75 2 — \6sm 2 + lOsmw 2 + I6(m 2 — mw 2 ) 2 ^jz 

+ s(3s 2 + 6sm w 2 -\4sm t 2 + Sm t 2 (m t 2 -m w 2 )y 2 + s 3 z 3 Df(z)J 
+ (z--z), (C2) 



. 2 12-N 2 +Af 2 z 2 m f 2 1 

27ia.a TTTTT^ — 



s N(N 2 -1) 6vr m w 2 l+z 



sz(l+z) +2(2(1 + z ) +z 2 ) [A (m 2 ) -A (m w 2 : 

- ( 1 + z) (is + (s + lm 2 - 2m w 2 )z s j B Q (s, m 2 , m 2 ) 
+ (2 + z) (s(l+z) - 2(m t 2 - m w 2 j)B~o(t, m t 2 , m w 2 ) 

- (s 2 + 4sm w 2 + S(m 2 -m w 2 ) 2 

1 — Z V 

+ 2s(s — 2m 2 + 2m w 2 )z + s 2 z 2 ^ 

x (Co(s,m 2 ,m 2 ,m 2 ) - (1 + z)Co(t , m 2 , m 2 , m w 2 ) j 

- -r~ \s(s-A(m t 2 -m w 2 )) 
1 -z L 

+ 2(s 2 — sm 2 + 2smw 2 + 3(m 2 — mw 2 ) 2 ^jz 

+ s 2 z 2 + 2(sm 2 + (m 2 -m w 2 ) 2 )z 3 Co(s,m 2 ,m w 2 \m 2 ) 

1 r3,/r2 2 o 2 2 

- — r s + os mw — 2s m t 

2(1 -z)L 

r\ r\ r\ sy 

+ \6smw (m\Y —m t ) + l6(mw — m t ) 

+ s(3s 2 + \2sm\y 2 — Ssm 2 + \6(m 2 — m\y 2 ) 2 ^jz 

(3s(s - 2m t 2 + 2m w 2 ) + Sm t 2 (m t 2 - m w 2 )}z 2 + s 3 z 3 ] D w (z) j 



+ s 
+ (z 



-z). 



(C.3) 



32 



References 

[1] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 327, 49 (1989) [Erratum-ibid. 
B 335, 260 (1990)]. 

[2] W. Beenakker, H. Kuijf, W. L. van Neerven and J. Smith, Phys. Rev. D 40 (1989) 
54. 

[3] R. K. Ellis and J. C. Sexton, Nucl. Phys. B 269 (1986) 445. 

[4] W. T. Giele, E. W. N. Glover and D. A. Kosower, Nucl. Phys. B 403, 633 (1993) 
[arXiv:hep-ph/9302225]. 

[5] J. H. Kiihn, A. Kulesza, S. Pozzorini and M. Schulze, Nucl. Phys. B 727, 368 
(2005) [arXiv:hep-ph/0507178]. 

[6] J. H. Kiihn, A. Kulesza, S. Pozzorini and M. Schulze, JHEP 0603, 059 (2006) 
[arXiv:hep-ph/0508253]. 

[7] J. H. Kiihn, A. Kulesza, S. Pozzorini and M. Schulze, Nucl. Phys. B 797, 27 
(2008) [arXiv:0708.0476 [hep-ph]]. 

[8] S. Moretti, M. R. Nolten and D. A. Ross, Phys. Lett. B 639, 513 (2006) [Erratum- 
ibid. B 660, 607 (2008)] [arXiv:hep-ph/0603083]. 

[9] S. Moretti, M. R. Nolten and D. A. Ross, Nucl. Phys. B 759, 50 (2006) 
[arXiv:hep-ph/0606201]. 

[10] W. Beenakker, A. Denner, W. Hollik, R. Mertig, T. Sack and D. Wackeroth, Nucl. 
Phys. B 411 (1994) 343. 

[11] J. H. Kiihn, A. Scharf and P. Uwer, Eur. Phys. J. C 45 (2006) 139 [arXiv:hep- 
ph/0508092]. 

[12] J. H. Kiihn, A. Scharf and P. Uwer, Eur. Phys. J. C 51 (2007) 37 [arXiv:hep- 
ph/0610335]. 

[13] W. Bernreuther, M. Fuecker and Z. G. Si, Phys. Rev. D 74, 113005 (2006) 
[arXiv:hep-ph/06 10334]. 

[14] J. H. Kuhn and A. A. Penin, arXiv:hep-ph/9906545. 

[15] J. H. Kuhn, A. A. Penin and V. A. Smirnov, Eur. Phys. J. C 17, 97 (2000) 
[arXiv:hep-ph/9912503]. 

[16] J. H. Kuhn, S. Moch, A. A. Penin and V. A. Smirnov, Nucl. Phys. B 616, 286 
(2001) [Erratum-ibid. B 648, 455 (2003)] [arXiv:hep-ph/0106298]. 



33 



[17] B. Feucht, J. H. Kuhn, A. A. Penin and V. A. Smirnov, Phys. Rev. Lett. 93, 
101802 (2004) [arXiv:hep-ph/0404082]. 

[18] B. Jantzen, J. H. Kuhn, A. A. Penin and V. A. Smirnov, Nucl. Phys. B 731, 188 
(2005) [Erratum-ibid. B 752, 327 (2006)] [arXiv:hep-ph/0509157]. 

[19] M. Beccaria, G. Montagna, F. Piccinini, F. M. Renard and C. Verzegnassi, Phys. 
Rev. D 58 (1998) 093014 [arXiv:hep-ph/9805250]. 

[20] P. Ciafaloni and D. Comelli, Phys. Lett. B 446 (1999) 278 [arXiv:hep- 
ph/9809321]. 

[21] V. S. Fadin, L. N. Lipatov, A. D. Martin and M. Melles, processes," Phys. Rev. D 
61 (2000) 094002 [arXiv:hep-ph/9910338]. 

[22] M. Beccaria et al, Phys. Rev. D 61 (2000) 01 1301; D 61 (2000) 073005. 

[23] A. Denner and S. Pozzorini, Eur. Phys. J. C 18 (2001) 461 [arXiv:hep- 
ph/0010201]. A. Denner and S. Pozzorini, Eur. Phys. J. C 21 (2001) 63 
[arXiv:hep-ph/0104127]. 

[24] E. Maina, S. Moretti, M. R. Nolten and D. A. Ross, Phys. Lett. B 570 (2003) 205 
[arXiv:hep-ph/0307021]. 

[25] G. Passarino and M. J. G. Veltman, Nucl. Phys. B 160 (1979) 151. 

[26] G. J. van Oldenborgh, Comput. Phys. Commun. 66 (1991) 1. 

[27] A. Denner, Fortsch. Phys. 41 (1993) 307 [arXiv:0709.1075 [hep-ph]]. 

[28] J. H. Kuhn and R. G. Stuart, Phys. Lett. B 200, 360 (1988). 

[29] S. Jadach, J. H. Kuhn, R. G. Stuart and Z. Was, Z. Phys. C 38, 609 (1988) 
[Erratum-ibid. C 45, 528 (1990)]. 

[30] S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 
510 (1998) 503] [arXiv:hep-ph/9605323]. 



34 



